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Abstract. 

We consider propagation of optical pulses under the interplay of dispersion and Kerr 
non-linearity in optical fibres with impurities distributed at random uniformly on the 
fibre. By using a model based on the non-linear Schrodinger equation we clarify how 
such inhomogeneities affect different aspects such as the number of solitons present 
and the intensity of the signal. We also obtain the mean distance for the signal to 
dissipate to a given level. 
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On the effect of random inhomogeneities in Kerr-media modelled by NLS equation 2 



1. Introduction 

In this article we consider the evolution of a complex electric field u(x,t) in a non- 
linear Kerr media which has constant dispersion and losses and, in addition, impurities 
at certain points x n , x n < x n+ i, which occur randomly on the fibre. We suppose that 
these loss elements cause the "input" signal u(x~,t) to abruptly decrease to an "output" 
value u(x+,t) = e~ 7n w(x~, t), where e~ 7n < 1 measures the dimming ratio and u(x~,t), 
say, denotes the limit value from the left. Assuming the validity of the self-focusing 
non-linear Schrodinger (NLS) equation as a model of ideal transmission [I] we find that 
the above situation must be described by a perturbed NLS equation which written in 
dimensionless units reads 



where the Dirac-delta terms account precisely for the amplitude decrease at impurities; 
further T > is the normalized loss coefficient. For the sake of avoiding extra 
mathematical difficulties we do not consider a compensated loss mechanism; this will 
be the subject of a future publication. We also remark that with minor changes our 
results may be applicable to other physically interesting systems such as Bose-Einstein 
condensates. 

It appears that while the effect of continuous random noise — or white noise — 
on NLS solitons has been well studied in the literature (see [21 El HI El E]) far less 
is known as regards the effects of sudden, discrete random perturbations. We intend 
to clarify how these inhomogeneities — which may be relevant for long-distance fibre- 
optic communication systems — affect the evolution of the pulse. We remark that 
perturbations involving delta masses also appear related to erbium-doped amplifiers and 
dispersion management, see 0E1E]- I n sucn a context, the positions of the amplifiers x n 
are deterministic and periodically disposed, while the strengths are constant 

and negative, 7„ = — Txi. Kodama and Hasegawa [10] generalize the latter ideas to 
a random context but, unlike us, maintain the amplifier interpretation and consider 
the distribution of the "intensity" of the signal only in the limit when both A n (here 
A n = x n — x n ^i > is the distance between impurities) and 7 n tend to zero. Thus while 
these ideas have some bearing with our work both the physical interpretation and the 
mathematical model are quite different. 

We will start our analysis of equation (pP) by considering that there are no 
deterministic losses, T = 0, since this case is simpler from a mathematical viewpoint: 
We show that upon performing a change of dependent variable the resulting formula 
can be piecewise related to the unperturbed NLS equation. Let us recall here that the 
classical NLS equation 



was first derived by Zakharov [TT] as an equation of slowly varying wave packets of 
small amplitude. He showed that despite its non-linear character the corresponding 




(1) 



n 



iQ x + Q tt + 2|9| 2 9 = 0, 9(0, t) = <p(t), 



(2) 
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initial value problem (IVP) can be reduced to a linear problem (the Zakharov-Shabat 
spectral problem) by the so called inverse scattering transform (1ST) — see [T2"] [T5] 
for general background on NLS equation and the 1ST method. Its interest has been 
further underlined by the realization that it also models the evolution of the complex 
amplitude of an optical pulse in a non-linear fibre [TJ. Applications of NLS equation to 
optical communications and photonics are nowadays standard [H [TJ [HI [15]. We devote 
section [2] to the study of the non-linear dynamics of the classical solitary waves within 
this regime, and we will show how impurities result in the appearance of radiation and 
general broadening of the signal. In particular, we find that solitons may be destroyed 
by the action of just one impurity. 

When r ^ equation ([T]) is no longer solvable in analytic way by 1ST; however 
we find — see section [3] — that the evolution of intensity, momentum and position of 
the pulse can be described precisely and that, under certain natural assumptions, 
their average values decrease exponentially due to the "impurities": concretely, we 
suppose that positions and strengths of impurities are statistically independent between 
themselves; we also suppose that in any interval [0, x] impurities are uniformly 
distributed (provided its number is given). Nevertheless the frequency and position 
of the pulse are not affected. 

In section H] we study the mean distance for the signal's intensity to attenuate to 
a given level due to the impurities. In applications, this level could be a recommended 
threshold value for reliability, say. To this end we formulate a linear integral equation 
that this distance satisfies and, by means of a Laplace transform, solve it. Results are 
discussed. 



2. Method of solution and the loss-less case 



Here we solve (jTJ given arbitrary sequences x n and 7 n with < x n < x n+ i and 7„ > 0. 
We perform the change of variable u(x,t) = ((x)v(x,t) where we require that ((x) 
depends only on space and has jump discontinuities at x n and that v(x, t) be continuous. 
By substitution we find that these functions must solve the equations 

iv x + v a + ( 2 \v\ 2 v = 0, (3) 
^ + rc(x) + £ (i - ^ n )S(x - x n )C(x-) = 0. (4) 

n 

It follows that ((x) is continuous on the intervals (x n ,x n+ i) wherein it solves equation 
01]) with no delta terms; further, it has jump discontinuities at the random points x = x n 
at which C( x n) = e ~ 7n (( x n)- Hence if N(x) is the number of defects on [0,x] we have 
that 

N(x) 

((x) = e~ s[x \ where S(x) = s + Tx + 7r ( 5 ) 

Alternatively, S(x) = s + Tx + YTj=ilji if x n < x < x n+i . (By contrast, for erbium- 
doped fibre amplifiers S(x) = —Tx — nTx\ if nx\ < x < (n+ l)xi, see |8j.) Thus S(x) is 
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a piece-wise linear function with initial value s and jumps at the random points x = x n , 
i.e., a pure random point process with drift, well known in the physics literature. For 
convenience we take s = hereafter, and until section HI There we will need to consider 
a more general situation where the starting value 5(0) is free. In figure [TJ we plot a 
sample of both S(x) and ((x) for a particular choice of the parameter set under this 
assumption. 




Figure 1. A sample path of ((x) (in red) and S(x) (green line) versus distance (in 
Km) showing the distance x for the energy to dissipate to half its initial value. We take 
a fibre with mean impurities distance < A„ >= A -1 = 1 km, a loss rate 0.02dB/Km 
and dispersion distance 50 km, i.e., T = 0.1 — which accounts for the seemingly linear 
behaviour between jumps. 

We shall now focus our attention in equation ([3]). We first consider the simpler case 
when the loss vanishes: T = 0. It turns out that, even though the resulting equation has 
random discontinuous coefficients, it can be piecewise reduced to an integrable equation 
whereupon we show how to obtain the evolution of an initial pulse (see [16] for related 
considerations). The reasoning in the rest of this section is essentially independent 
of the sequences A n = x n — x n _i and j n . Nevertheless we shall suppose that both 
are sequences of positive, independent, equally distributed random variables and that 
A n and 7 m are also independent for all n, m. Note that all these assumptions are 
physically well founded as they imply, say, that the knowledge of the position of a given 
impurity does not provide any information on the location of the remaining ones. The 
further assumption that A n is exponentially distributed: Pr ^A n > xj = e~ Xx where 
A =< A n > _1 is a certain parameter, is natural from physical principles. It has several 
fruitful consequences as then there follows that the number N(x) of impurities that 
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occur on [0, x] has Poisson distribution with parameter Ax and that they are uniformly 
distributed on the interval. It further implies the memory-less property: the distribution 
of impurities on (x, x + Ax] remains unaffected given that none was observed on [0, x]. 
By contrast, we consider here a general probability density function (PDF) h(y) of 7 n : 
Pr (y < In < V + dy) = h(y)dy. 

For the sake of being specific let us consider the case when the initial data is just 
a solitary wave pulse: v (0,t) = 2^sech (2i]t) e 2 *^' = ^^(t) where the real parameters n 
and £ give, up to a constant, the wave's amplitude and the carrier velocity. Note that 
up to the first impurity v^°'(x, t) = v(x, t), < x < xi, solves the I VP 

iv<® + 4 0) + 2|t/°)|V°) = 0, u(°)(0,t) = <p {0) (t). (6) 

This is the standard IVP for NLS equation and hence the solution for < x < x\ is the 
classical soliton 

u(x,t) = i/ 0) (x,t) = 2r/sech(27 ? (t-4ex))e i [ 2 «' +4 ^ 2 -« 2 H. (7) 

As commented, we continue this solution to the interval Xi < x < x<i by requiring v(x, t) 
to be continuous at x = x\. This requirement fixes v^(x,t) = v(x,t), x\ < x < X2, to 
satisfy the non-linear partial differential equation 

iv W +V W + 2e- 2 ^\v^\ 2 v^ = 0, with 

v {1) (x u t) = 2r] sech (2 V (t - 4^xi)) e 4 2 €'+ 4 (" 2 ^ 2 )^]. 

Remarkably this equation can be reduced again to NLS: by using both temporal and 
translational invariance of NLS equation one can prove that 

e 7l w(x, t) = v w (x, t) = e 4i(v2+e)xi+ ^e(x -x u t- 4£xi), 

where 6(-,-) is the solution to the NLS equation ([2]), with data at x = given by 
6(0, t) = e" 7 V (0) W- Notice that, unlike v(x, t), u(x,t) is not continuous at x = x\. 

The determination of the specific form of the solution requires solving a linear 
spectral problem. The procedure is awkward but fortunately the solution's main features 
may to a large extent be determined avoiding these complexities. We note that if 71 7^ 
the solution that evolves from data O(0, t) = e~ 7l <^^°- ) (t) is no longer a soliton but a 
complicated pulse that may contain radiation, in addition to the soliton. The former 
component has a much weaker rate of decay than the later; concretely, it decays as the 
corresponding solution for the linearized Schrodinger equation (i.e. as t~ 1/2 , see [12]). 
Further, if 7x > 1.41 the arriving soliton at x = Xi — cf. equation ([7]) — is destroyed 
by the action of the fist impurity after X\. Hence the resulting configuration for 
x > xi consists solely of radiation. To be specific, suppose that the jump PDF h(-) 
has exponential distribution with mean 1/er. Then, after the first impurity the soliton 
disappears with probability bounded below by Pr(7i > y) = e~~ ay , where y = 1.41. 

f We adopt the convention and terminology of standard NLS theory wherein t is space and x a temporal 
variable, a situation opposite to that that occurs in Optics. 

I This stems from the fact that the condition S 2 io (2S) < 1 on the initial data guarantees that no 
solitons will be formed upon evolution [12l [17] . Here Jo(-) denotes the modified Bessel function of zero 
order and 5 = \@(0,t)\dt = Tre" 71 . 
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Finally, we mention that by using similar ideas one can extend the solution to 
x > x n by solving (J3J) with data ,t), where as before v^ n \x,t) denotes the 

general solution v(x,t) restricted on x n -\ < x < x n . Translation invariance allows 
one to reduce this to NLS equation with new data which involves a contraction factor 
e -(7i+-+7n)_ Eventually, this dimming of the initial signal results in a disappearance 
of the starting solitons into radiation, an indication that, as a result of impurities, 
broadening of the signal takes place. We skip the mathematical details. 

3. General case with deterministic loss and impurities 

When T > equation ([3]) can be mapped into the so called dispersion- managed 
NLS equation, which, unfortunately, is not solvable in analytic way, neither by using 
1ST nor by any other method. It is then remarkable that the evolution of the main 
physically observable functional can be discerned in an exact way. Consider the 
following quantities 

/oo 
\u(x,t)\ 2 dt, 
-oo 

/oo 
u(x,t)u t (x,t) dt, and 
-oo 

/oo 
t\u(x,t)\ 2 dt, 
-oo 

where M(x) and P(x) are the (accumulated) intensity and momentum of the signal 
at a position x, while Q(x)/M(x) = T(x) is the pulse position. The functional 
P(x)/M(x) = Q(x) is interpreted as the pulse-centre frequency. The singular nature 
of the delta terms prevent us from determining the relevant evolution by manipulating 
equation ([T]). Nevertheless, one can rely again in the decomposition u(x,t) = ((t)v(x,t) 
and use equation ([3]). Then, proper manipulation of the latter expressions yields that 

M(x) = M(0)e- 2S{x \ 
P(x) = P(0)e- 2S(x \ and 
Q(x) = [Q(0)-2P(0)x}e- 2S ^. 

Thus the effect of the presence of impurities results in the addition of a multiplicative 
random factor e~ 2S<yX ^ in both intensity and momentum. Note however that Q(x) = fi(0) 
and T(x) = [T(0) — 2Q(0)x], and hence that inhomogeneities have no effect whatsoever 
on position and frequency, a fact that accords with the physical intuition. 

It is therefore of interest to evaluate the mean amplitude's decrease. We do so 
by first assuming that previously n defects have occurred: N(x) = n. Let E denote 
statistical averaging and E,(^ 2 (x)\N(x) = n^j be the mean value of ( 2 (x) knowing that 
exactly n jumps have occurred on [0, a;]. Note that given this information one has 
S(x) = Tx + Yl n j=\ If i- e -' om y the uncertainty regarding the value of the 7/s remains 
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but not that associated with the number of summands N(x). In view of the assumed 
statistical independence we have that the mean factorizes as 

E ( C 2 (x) | N(x) = n) = E i^e- 2Tx f[ e" 2 ^ j 

n 

= e- 2rx l[E(e- 2 ^ =e~ 2rx Ql 
i=i 

where Q r = E [exp ( — rjj)] = J °° e~ ry h(y)dy < 1 is the Laplace Transform of the 
jump-size PDF. The mean intensity is obtained by further averaging with respect to the 
number of impurities: 

E[M(x)} = M E[C 2 (x)] 

= M Q J2 , 6 E(( 2 (x)\N(x) = n) = M e- [2r+A(1 - Q2)1 -, (8) 

n=0 U ' 

where M(0) = M and we used that if Aj has exponential distribution, i.e., if 
Pr(A j >x)= e~ Xx for some A > 0, then N(x), the number of defects on [0, x], is Poisson 
distributed: Pr (N(x) = n) = (Xx) n e~ Xx jn\. Hence we obtain that the existence of 
defects implies an additional exponential decrease in the field's intensity and momentum 
at a rate 2A(1 — Q 2 ), an effect which might result in the degradation of the bit patterns. 



4. Mean half life 



A natural related problem of interest is determining the distance x at which M(x) 
dissipates from a starting value Mq to a given level M\, i.e., such that M(x) = Mi. For 
convenience we set Mi = M e~ 2b and hence require S'(x) = b. This distance could be 
considered as a threshold value below which the signal is no longer reliable (it gives the 
mean half life of the signal if M = 2Mi). In the deterministic case (A = 0) this distance 
follows inverting Mi = M exp ( — 2rx) as x = ^ log ^ . When inhomogeneities are 
present x is a random variable whose mean is not obtainedbj inverting equation (JSJ) — as 
it might have been naively thought. Instead, we reason as follows: call x s , see figured! 
the (random) distance that takes for the generalized process S(x) with initial value 
S(0) = s — cf. equation (jSJ) — to go beyond the level b. It turns out that X(s) = E(x s ) 
satisfies the linear integral equation 

X(a) = — |^ + fjy dle^ s+l ~ b) J dyX(y + b - l)h(y), (9) 

where g = ^jr- and we recall that h(x) is the density of 7„. Q 

* We sketch the derivation of this integral equation (see [18] for a similar derivation in a financial 
context). With S(0) = s there are three possibilities for the future evolution: If the first jump satisfies 
x\ > g then S reaches the level b at x = g. If this is not the case and if the jump at x\ satisfies 
s + Txi + 71 > b then the process goes past b at x = x\. Otherwise the process still remains within 
[0, b) at x = x\ and starts afresh with an initial value S(x\) = s + Tx\ + 71 < b (hence the process will 
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This equation can be solved in a closed form by Laplace transformation. We 
consider again the case corresponding to a jump PDF also exponential with mean 
a -1 =< 7 n >, i.e., h(x) = ae~ ax where a > 0. If k = A + err, Laplace transformation 
yields the solution to (P) as 



aTg A 



K 



-Al-e 



—kq 



The mean distance for the amplitude to decrease to Mi follows letting s 
b- 1 



2 log f| as 



E <*™-2(FTAM log f 4 



1 



/Mi 
I M 



and 



(10) 



If A 



we recover the deterministic limit above: x = 2f log Iff. Note how the 
incorporation of impurities corrects this formula in a significant way, cf. equation ffTOjh 
Another interesting limit is that of vanishing deterministic loss rate, T = 0. The mean 
attenuation distance can only be accounted to the presence of impurities and reads 
E(x) = j + ^logffj- The first term is the mean time for the first jump at X\ to 
happen; the logarithmic correction corresponds to the mean time to go beyond the level 
b after the first jump. Actually, this rate rules the mean dissipation distance whenever 
M >> Mi and A/a >> T. In figure [2] we perform a plot of this function. Note how, 
by contrast, the distance implied inverting equation (J8]), namely 

A 



E(M(x)) 



2 r + 



M exp 
A 



-2x r + 



a 



and therefore 



a + 2 



~i -i 



log 



Mo 
Mi' 



(12) 



deviates from the correct result, equation f lTU]) . and fails to capture the sharp behaviour 
occurring for Mi Mq. The error increases as T decreases. 



5. Conclusions 



We have analyzed how the existence of randomly distributed impurities affects the 
evolution of an optical pulse in a non-linear Kerr media with constant dispersion and 
loss. We suppose that the unperturbed situation is described by NLS equation. When 
the deterministic loss vanishes it is shown by changing the dependent variable that 
the resulting equation can be piecewise related to the unperturbed NLS equation. The 
effect of impurities in the non-linear propagation is pinpointed. In particular we address 
the issue of how they affect the initial solitons and the possibility to dissipate them 

exit [0, b) at x\ + x s+r:ci+71 ). Upon appropriate rearrangement this reasoning leads to 
x s = g0(xi -g)+x 1 9(g-x 1 )+x s+rxi+ ' 11 9(b-s-Tx 1 -' ll )d(g-x 1 ). 
Averaging this relationship yields with further manipulations equation 
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2.0 




1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 

M /Mj 

Figure 2. Mean distance in terms of Mq/Mi for A = 2.0, a = 3.0 while T = 0.5 (red 
line) and T — 0.05 (blue one) as follows from (fTU)) . Note how in the latter case X(0) 
jumps an amount < A„ >= 0.5 right after the origin. The green and magenta curves 
are the (incorrect) mean distances implied by equation (|12[) with the above parameters. 

into radiation. In the general, non-solvable r ^ case we show that while impurities 
do not influence the frequency and position of the signal they induce an exponential 
decrease of the main physical observables intensity and momentum and hence a general 
degradation. We also determine the mean half life or mean distance for the signal to 
dissipate to a given threshold value. We find that this distance satisfies a certain integral 
equation. Its analysis shows that impurities result in an important decrease in the mean 
dissipation distance. To overcome these effects the addition of amplifiers is in order. The 
introduction of such a device and the relevant statistical implications will be considered 
in a future publication. 
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